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ABSTRACT 

> 

The D = matrix model is reformulated as a nonlocal quantum field the- 
^ . ory in two dimensions, in which the interactions occur on the one-dimensional 

^ ■ line of hermitian matrix eigenvalues. The field can be thought of as a fluctuation 
in the potential V, and is conjugate to the density of matrix eigenvalues which 
appears in the Jevicki collective field theory. The classical solution of the field 
equation is either unique or labeled by a discrete index. Such a solution corre- 
sponds to the Dyson sea modified by an entropy term. The modification smoothes 

• , 

the sea edges, and interpolates between different eigenvalue bands for multiple-well 
potentials. Our classical eigenvalue density contains nonplanar effects, and satis- 
fies a local nonlinear Schrodinger equation with similarities to the Marinari-Parisi 
D = 1 reformulation. The quantum fluctuations about a classical solution are 
computable, and the IR and UV divergences are manifestly removed to all orders. 
The quantum corrections greatly simplify in the double scaling limit, and include 
both string-perturbative and nonperturbative effects. The latter are unambiguous 
for V bounded from below, and can be compared with the various nonperturbative 
definitions of these theories proposed in the literature. 
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1. Introduction. 



The study of discretized two-dimensional quantum gravity (matrix models) has 
received much attention in recent years, thanks to the uncovering of the solvability 

[1 — 8] 

of a range of such models in the double scaling limit . A matrix model describes 
the statistical mechanics of random surfaces in the large- N limit, where N is the 
dimension of the matrix (we restrict our attention to hermitian matrix models). 
The procedure for solving these models consists of three conceptual steps. In the 
first step, integrating out the angular modes of the random matrix leaves us with its 
N real eigenvalues as the dynamical variables. Then, the large-N limit is performed 
via a WKB-like procedure^. And finally, as N is taken to infinity, the matrix-model 
potential must approach one of its critical values in accordance with the double- 
scaling limit. The averages of various quantities over random surfaces can then be 
found as functions of the string coupling, which is held fixed in this limit. To any 
given order in the string-perturbative expansion, the matrix model results may be 
compared with results from the corresponding continuum models, which consist 
of the quantum Liouville theory coupled to various matter fields.' 9,101111 141 . Such 
comparisons have been carried out for low genus [1 ° 181 . 

In addition to perturbative series in the string coupling, one finds in these 
models also nonperturbative effects' 191 . In the D = 1 model, the leading such 
effect comes from the tunneling of a single matrix eigenvalue (=fermion ) out of 
the potential well, or between different potential wells. The contribution to, say, 
the free energy will behave as exp(— const / g string) f° r l° w string couplings. In 
ref.19 it was pointed out that for the D = models, as well, there are saddle-point 
configurations where a single eigenvalue leaves the potential well"'", and that these 
configurations give rise to nonperturbative effects of the same form. However, for 
even-A; multicritical D = models, including pure gravity itself, there are ambigu- 

f For the D = 1 model, this is a true quantum-mechanical WKB procedure for N independent 

fermions in a potential well. 
| In the D = case, this is the effective potential, namely the matrix-model potential plus 

the mean potential due to the Coulomb repulsion of the other N-l particles. 



2 



ities in denning the theory nonperturbatively, which can be traced to the fact that 
the critical quartic potential is unbounded from below. Various nonperturbative 
definitions have been proposed in the literature' 20 311 . For the even-k multicritical 
models, these definitions can and do disagree with one another, since the per- 
turbative series is not Borel-summable. But even for the well-defined models, a 
systematic derivation of the nonperturbative physics directly in D = has so far 
been lacking. 

A further line of development has been the reformulation of the D — 1 matrix 
model as a string field theory. In matrix language, the correlators that one calcu- 
lates in this model are of any number of boundary operators, each such boundary 
having an arbitrary length. Each boundary is in a single slice of the embedding 
dimension . The same information is contained in the n-point functions of the den- 
sity of matrix eigenvalues; it is essentially this density which has been suggested as 

[32] 

the string field . It is a function of A, the matrix eigenvalue, and the embedding 
dimension; hence the field theory is two-dimensional. In ref. [32], Das and Jevicki 
used collective coordinate techniques to transform the quantum mechanics of N 
particles (in a bosonic formulation) into two-dimensional quantum field theory. 
The kinetic term in the action is that of a massless field, which corresponds to the 
massless tachyon in the continuum effective field theory [12]. Some aspects of this 
correspondence remain unclear — for instance, the identification of matrix eigen- 
value with (a function of the) Liouville zero-mode. Nevertheless, much progress 
has been made in understanding the Das- Jevicki collective field theory, making 
string-perturbative calculations with it, and comparing the results with those ob- 
tained via other methods 133 ' 341 . Nonperturbative effects, appearing in the form of 
solitons and instantons of the D — 1 field theory, have also been investigated 135 ' 361 . 

Another approach to D = 1 field theory has been to retain its formulation as N- 
fermion quantum mechanics, but second- quantize the fermions and then bosonize 
them 137 391 . This approach is equivalent to the Das- Jevicki formulation. 

§ Otherwise, the angular matrix variables are excited and the calculations cannot be done 
using current methods. 
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Although the collective-field method has also been applied to the D = matrix 
models' 40 431 , it seems to us less compelling there than for D — 1, because of its 
unusual kinetic term, and although a perturbative scheme has been outlined [44,4j1 , 
this program has not, to our knowledge, been carried out as far as the corresponding 
D = 1 field theory 1 . 

In this paper, we develop an alternative field theory formalism for the D = 
matrix models. Our field theory is two-dimensional, but the extra dimension is 
an auxiliary one and drops out of the formalism at some stage. What is left is a 
nonlocal quantum mechanics, with a function of A playing the role of time. This 
is interesting in view of the idea [32] that the matrix eigenvalue is related to the 
Liouville zero-mode. The Jevicki-Sakita D = collective field theory is also a 
nonlocal quantum mechanics. We have found, however, that our formalism has 
the following desirable properties, which make it worth pursuing: 

A. The entropy term* in the action, which in the collective field approach 
appears through the Jacobian, is for the first time incorporated into the classical 
solution. This results in smoothing of the Dyson sea edges, interpolation between 
different eigenvalue bands for multiple-well potentials (determining their relative 
population), and an unambiguous treatment of nonperturbative (instanton) effects, 
directly in a D = field-theory framework, for multicritical models that can occur 
for a bounded-from-below potential. 

B. Our classical eigenvalue distribution, p(A), satisfies a local nonlinear 
Schrodinger equation, with 1/N playing the role of U. The Schrodinger potential, 
Vi(A), appearing in this equation is similar to, but distinct from, the Marinari- 
Parisi D = 1 reformulation. On the other hand, in the planar limit it has an exact 
corresposondence with the Dyson sea effective potential, mentioned above. 

% In ref. [45] , Lcchtcnfcld uses our equations to set up a perturbative scheme. But since he 
performs the genus expansion at an early stage, one is left with manifest divergences, and 
in addition nonperturbative information is lost. We avoid these problems, as will be seen 
below. 

* This term is N Jplnp dX, with p(X) the eigenvalue density. 
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C. Not only are quantum corrections computable, but most of them vanish in 
the double scaling limit. The only quantum corrections which survive this limit, 
apart from the semiclassical functional determinant, are a set of Feynamn graphs 
which can be exactly summed. 

The rest of this paper is organized as follows. In section 2, the matrix-model 
partition sum is recast as that of a massless two-dimensional field, A(r), with (non- 
local) interactions confined to an infinite line (the 'eigenvalue axis'). The equation 
of motion is written, and its classical solution is expressed as an eigenvalue distri- 
bution, p(A), and shown to be either unique or labeled by a discrete index. The 
two-dimensional equation of motion becomes a one-dimensional integro-differential 
equation on the eigenvalue axis; this in turn leads to a weaker** local nonlinear 
Schrodinger equation. The interpretation of this equation is discussed, as well as 
the nature of its solution and how it is determined uniquely (up to the possible 
discrete index). 

In section 3, the quantum fluctuations about a classical solution are studied, 
and an expression derived for the partition sum in terms of the propagator of the 
A field (the 'conjugate field'). The IR and UV singularities are seen to cancel in a 
trivial way, to all orders. The one-dimensional Green's equation for the propagator 
is presented, as well as a partial solution. Section 4 is devoted to a discussion of 
the double scaling limit, including nonperturbative effects***. We see there that 
the quartic k = 2 classical solution is unique at the level of string perturbation 

** 'Weaker' because the Schrodinger equation has more solutions than the integro-differential 
equation. 

*** Most of the discussion in section 4 is limited to the case of the quartic potential as it 
approaches its k — 2 critical point. Since neither the theory nor our formalism are well 
defined for this potential, only the string-perturbative series discussed in that section is 
meaningful. Therefore, our description of an instanton ansatz for this model, described in 
sec. 4, is not necessarily more reliable than any of the previous nonperturbative definitions 
of pure gravity referenced above. Nevertheless, it is reassuring that our exponentially- 
suppressed tunneling factor, cq.(39a), agrees with that appearing in other approaches (it is 
the prefactor that is ambiguous). The main point of our instanton calculation is, however, 
that the conjugate-field formalism applies to any well-defined realization of a critical model, 
and a similar instanton calculation for such a realization would be rigorous. 
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theory. In section 5 we restate our conclusions. Some mathematical details are 
reserved for the appendix. Many results, presented or stated without proof in this 
paper, will be exposed more fully in a follow-up publication, to appear soon. 

Finally, an explanatory note is in order. An early version of this preprint, 
containing most of the results (with the notable exception of point C above), has 
been circulating privately since November 1990. That earlier version has sometimes 
been confused in the literature with SLAC-PUB-5262, an altogether distinct work. 

2. Conjugate Field Formalism and Classical Solutions. 

The partition function of the D = matrix model is, 



Where V is the matrix-model potential, which we assume to be a polynomial. 
We may formally rewrite (C a divergent number) 



Which is the path integral over a massless field A in two dimensions, with iV 
point charges on the x = ('eigenvalue') axis, x is an auxiliary dimension, and 
A(X) = A(X0). We shall use r to denote a general point (X,x) in the plane. 

This path integral is beset by infrared and ultraviolet divergences. The UV 
divergence is regulated by smearing the point charges; A(Xj) is replaced by 



in eq. (2). The infrared divergence is regulated by introducing a uniformly charged 
circle in the plane, centered about the origin and with a large radius L. The total 
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A 3 +e 




(2a) 



charge of this circle is — N, which screens the N point-charges. All dependences 
on e and L will cancel in a simple way. In what follows, we will mostly ignore the 
need to regulate the path integral; a more careful treatment reveals that this naive 
approach is the correct one. Introducing the normalized charge density, 



pW = ^J2 6 ( x ~ x ^ ( 3 ) 

i 

We may instead think of Z as a path integral over the overcomplete variables 
{p(A)}, with integrand 



exp(-iV 2 J d\p(\)V(\) + iV 2 J J d\dnp{\)p{n) In |A - ^|) 



(4) 



This is the collective-field approach of Jevicki and Sakita. But we prefer to 
use the conjugate field A(r) as our dynamical field, since it has a standard kinetic 
term and trivial Jacobian. The other merits of the conjugate field theory have 
been listed in points A, B and C of the introduction. 

Combining eqs. (1) and (2), and using (2a) to regulate the UV divergence in 
C, we obtain the following expression for the partition function*: 

Z{V} = C l (L,N)e~ N J[dA](l{A}) N exp(-^=^dsA-^ j \dA) 2 d 2 r) (5) 



I{A} = J d\exp(-NV{\) + iV^A(X)) (6) 

In equation (5), the linear term is an integral over the charged circle, and C\ 
is IR divergent, UV finite, and V-independent. All integrals over A, here and 
below, range over the whole real axis, unless the limits of integration are explicitly 
indicated. 

* We define the measure [dA] to include the infinite, but {L, e, N, V(X)} independent, factor 
of (det'd 2 ) 1 ^. 
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The euclidean action is, apart from the IR-regulator source term, 

S{A} = ^J {dA) 2 d 2 r - N\nI{A}. (6a) 

This is a nonlocal action, with the interactions occuring on the eigenvalue axis. 
Away from the axis, A is a free massless field. 

i\fA^A/N can be thought of as a quantum fluctuation in the matrix-model 
potential. In addition, A is conjugate to p in the thermodynamical sense. 

The classical equation of motion for the field A is^ 

d 2 A(r) = -iV^N Pl (X)5(x) (7a) 

where 

pi(A) = j exp(-W(A) + iV^A(X)) (7b) 

Note that although the action is nonlocal, the equation of motion is local, but 
with a new coupling /, to be determined self-consistently. 

The boundary condition for the conjugate field A at infinity is, 

2N 

A(r) ~ —i In jr| + const (7c) 

V47T 

since that is the electrostatic potential far from the N point-sources (eq.(2)). This 
boundary condition allowed us to freely integrate by parts the variation of the 
kinetic term in the action. 

f We omit a step here, related to the infrared regulator: the field A should be shifted by a 
constant, which is the constant potential due to the charged circle in its interior. 
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We denote a solution of eqs.(7) by 



classical V 



iA s (r) 



{Id) 



and define the classical charge (eigenvalue) density on the eigenvalue axis as p at 
the classical solution: 



Where I s is the value of I{A} at the classical solution. 

Henceforth we shall adopt the definition eq.(8) for p, instead of eq.(3). From 
here on we shall work with p instead of A s ; They contain the same information, 
although p is defined only on the eigenvalue axis"'" . 

Let us assume that the potential V has been chosen to be bounded from below. 
In that will be shown, p(A) is unique or labeled by a discrete index. 

In section 4, we shall see that quantizing A{r) about the classical configuration 
determined by p, gives rise to a unique string-perturbative expansion for physical 
quantities in the double scaling limit (abreviated henceforth as d.s.fi ), provided p 
has support in a single Dyson sea (band) in the planar limit. Thus, the possible 
discrete non-uniqueness of p will be revealed either at the nonperturbative level, 
or for multiband solutions. However, any continuous non-uniqueness of multiband 
solutions, which results from the freedom to adjust the relative populations of 
different bands, is eliminated by the nonperturbative effects ('tunneling'), as we 
shall see. 

| Just as, in the analog 2d electrostatic problem, the charge distribution on a plate together 
with the boundary condition at infinity, determine the potential throughout an otherwise 
empty space. 

§ In order to properly analyze the d.s.l., a multiple- well potential, bounded from below, 
must be used. In this paper, when we specialize to a particular V(X), it is the quartic 
A 2 /2 + <?A 4 , which is not bounded from below at criticality. This is not important for string 
perturbation theory; work on nonperturbative effects in which a multiple-well potential is 
used, is currently in progress. 




s 
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Both the equation of motion and boundary condition for A s (r) are real, and 
thus each p(A) in the discrete set of solutions is either real, or the complex conju- 
gate of another solution. The possible imaginary parts of solutions p can only be 
revealed nonperturbatively; in any case, all physical quantities are real. 

Evaluating the action at the classical solution, we find that the L-dependent 
prefactor in eq. (5) gets canceled: 



ldassical{V} =const e N exp(-iV 2 / dXpV 

(5a) 



+ N 2 J J d\dpp{\)p{p) \n\X-p\-N J d\p\np) 



The constant depends only on N, and is hence irrelevant. 

As we discuss below, the UV divergence in eq (5a) is cancelled by the semiclas- 
sical (determinant) factor, whereas the higher quantum corrections to Z are both 
IR- and UV-Gmte. 

Notice the O(N) correction to the free energy that occurs already at the classi- 
cal level. This correction has a simple interpretation: It is a combinatorical factor, 
the entropy of the classical solution (see refs. [40], [41], [42], [43]) . Of course, the 
sum of all O(N) corrections to the free energy must vanish, since only even powers 
of N appear in the topological genus expansion' 461 . 

p satisfies the integral equation, 



\/(A) + ^ln(p(A)/ s )=2 J dp\u\\-p\p{p) (9a) 

which follows from eqs.(7), in conjunction with the fact that the 2d free Green's 
function is ^ In \r— r'\. Upon differentiation w.r.t. A, eq.(9a) becomes the following 
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one-dimensional integro-differential equation: 



V'(X) + ^^ = 2H(p(X)) (9b) 
Where H denotes the Hilbert transform: 

WW) - / ^ (9c) 

for any function f . 

In addition to satisfying the Id integro-differential equation, p must be nor- 
malized to unity (by eq. (8)): 



J d\p(X) 



(9d) 



The eq. (9b) reduces to the integral equation of ref.46 in the planar (N — > oo) limit; 
the extra term can be physically understood as due to the variation of the entropy 
term in the exponent of eq.(5a). Its effect is to smear the edges of the Dyson sea 
(or seas, for a multi-well potential V), so they become transition regions, rather 
than singularities as in the D = and D = 1 collective field theories of Jevicki et. 
al. We will further discuss the transition regions in section A m . 

For N >> 1, it is easy to see from eq.(9b) that outside the sea, and beyond 
the transition regions, p(A) is well-approximated by 

p(X) » const A^e-^W+^/A)] (9e) 

Thus, two conclusions can immediately be drawn. Firstly, a normalizable classical 
p exists only if V(A) is bounded from below; and secondly, p(X) vanishes nowhere, 
so that multiple eigenvalue bands ('seas') are related by a tunneling effect. This 
effect serves to determine the relative population of multiple seas, and is responsible 
in general for string-nonperturbative effects. 



* Such integrals are understood, here and below, to be principal- valued. Note that we use a 
nonstandard normalization in our definition of H. 
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We note that, as our classical solution depends on N and hence includes some 
higher-genus effects, the terms 'planar' and 'classical' are not synonymous in our 
formalism. This is also the case for the Das-Jevicki field theory. 

Once p is known, the corresponding I s is determined uniquely by eqs.(9); how- 
ever, I s drops out of the formalism from here on. 

In the appendix we use eqs.(9b-d), and properties of the Hilbert transform, to 
derive the following nonlinear Schrodinger equation: 

(-^^ + ^i(A) + 7 rV)(p- 1/2 )=0 (10a) 
where V\, henceforth called the 'Schrodinger potential', is the polynomial 

V 1 (\) = \(V) 2 + ±V» + P(\) (106) 
P is a polynomial whose coefficients are moments of the charge distribution: 



m ^j d , M m^m (10c) 

This Schrodinger equation is local, except for the coefficients of P{X), which are 
determined self-consistently. For example, in the case of quartic V, 

V(X)= 1 -X 2 + g\ 4 (lOd) 

the polynomial P(X) becomes: 

P(X) = -1 -4<?(m 2 + A 2 ) (lOe) 
where the eigenvalue moments are defines as 

m n = J X n p(X)dX (10/) 
and we have used the symmetry p(A) = p(— A), which follows from eq. (9b) and 
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the symmetry of V(A) . The nonlinear Schrodinger equation has several points 
of interest, which we now discuss. Firstly, in it 1/N plays the role of Planck's 
constant, and the tunneling effects implied by eq.(9e), can now be seen to be 
similar to quantum mechanical tunneling. This is demonstrated in section 4, using 
the WKB approximation. In the exterior of the Dyson sea(s), p is exponentially 
suppressed for large N, and (10a) approximates a linear Schrodinger equation, a 
strange feature of this correspondence is that the 'wavef unction' here is 1/y/p, the 
inverse of the intuitive yfp. Another feature, probably closely related, is that the 
V" term in our Schrodinger potential Vi, has the opposite sign compared with the 
D = 1 potential arising from the Marinari-Parisi reformulation [24] . 

The transition regions interpolate between the WKB solutions outside and 
inside a given Dyson sea, similar to the role of the Airy function in ordinary 
(linear) WKB (see section 4). In the sea interior, the p 2 term in (10a) is no longer 
negligible, and in fact there the leading WKB approximation is 

p(A) ~ — a/— Vi(A) (inside a sea) (10g) 

7T 

which is just the generalization of Wigner's semicircle law to arbitrary potential V 
[46] . 

In the planar limit and outside the Dyson sea, Vi(A) has a physical interpre- 
tation: 4Vi(A) is the square of the gradient of the effective potential, that is, the 
square of the total force exerted on a single eigenvalue, due to the potential V(A) 
and the repulsion of the other (N — 1) eigenvalues. 

f For simplicity, we assume throughout that V(X) = V(— A). Actually, this only implies a 
symmetric p if p(X) is unique, but we will ignore this complication here, since the entire 
analysis can be easily redone without the symmetry assumption, and besides the possible 
nonsymmetry will not affect perturbation theory. 

| The relation between the two facts is, that in order to get the Marinari-Parisi D = 1 
Schrodinger equation, one chooses as wavefunction ip{M} = cxp(— ^ tr V(M)), with M the 
original random matrix. Thus, if the Van der Monde were to be ignored, the one-particle 
wave-function would be the square root of p. If we choose the inverse wave function, the 
sign of V" in the Schrodinger potential agrees with ours, rather than with that of Marinari 
and Parisi. 
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As stated in the introduction, the Schrodinger equation is weaker than the 
integro-differential equation. This is because (10a) is a second-order differential 
equation, so for given V\ its solution p(A) has two free continuous real parameters. 
By eqs.(10), V\ has one additional unknown, m% but we must impose the two 
self-consistency conditions, for rri2 and for uiq = 1. This leaves us with a single 
undetermined real parameter . However, as is proven in the appendix, the integro- 
differential equation allows no zero modes that preserve the normalization condition 
(9d). Hence, the free real parameter can only assume discrete values. 



3. Quantum Corrections. 

We next address the quantum corrections to Z c i ass i ca i ( eq. (5a)), needed in 
order to regain the full partition sum Z in eq. (5). Let us separate the field A into 
its classical and quantum pieces, 

A(r) = -tA s (r) + A q (r) (11) 

and also separate out the quadratic part of the action: 

S{A} = S classical + S T {A q } + IJ J d 2 rd 2 r'A q (r)K(r, r')A q (r') (12) 

Here Sj is the interacting piece, consisting of the terms of order three and 
higher in A q . K is the inverse propagator of the quantum field in the background 
of the classical solution: 



§ Note that the usual WKB procedure for bound states, does not apply for the Dyson sea. 
In the usual procedure, one imposes that the wavefunction component which blows up 
exponentially at spatial infinity, vanishes. But here the wavefunction is ip = p^ 1 ^ 2 , so in 
fact p normalizability requires ip to blow up exponentially, and no useful information results 
from this condition. 
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K(r, r') = -(d r ) 2 S(r - r') + 4irN5(x)5(x'){p(X)5(X - X') - p(X)p(X')} (13) 



The only zero mode of K is the constant function^, so we fix that by defining 
our space of configurations A q to satisfy A q {r) — > as \r\ — > oo. This renders K 
nonsingular*. 

Now, Sj depends only on the value of A q on the x = axis, so we define the 
one-dimensional field 

q(X) = V^[A q (X0) - J p(p)dpA q (pO)] (14) 

and denote 

Si{q} ee 5 7 {AJ (15) 
We then have, by eqs.(6),(6a),(8),(ll),(12) and (14): 

^{g} = -A r ln(|p(A)rfAe^ A )) ivg (16) 

where NQ denotes the non-quadratic part. The full partition function, including 
quantum corrections, is given by 

Z{V} = j: Z classical{V}(^^y 1/2 (e- S ^) (17) 

where the expectation value is a Gaussian average, i.e. evaluated via Wick's the- 
orem, the sum is over the discrete set of classical solutions, the functional deter- 
minants are subject to the boundary condition rim| r |_ KX) A q (r) = 0, and for each 
classical solution, Z c i ass i ca i is given by eq.(5a). 



This is equivalent to the fact, proven in the appendix, that p(X) has no normalization- 
preserving zero modes. 

* This choice of boundary condition follows from the fact that both A(r) and —iA s (r) satisfy 
the boundary condition (7c). 
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The two-point function, to be used in the Wick expansion, is 



( q (\) q (\'))=4n(H(\,\') + ^ N ) (18) 

where H(X,X') is the restriction to the eigenvalue axis of H(r,r'); the latter is 
defined uniquely by the following properties — 



d?H(r, r') = 47rNp{\)5{x)H{r, r') + 5(r - r') (19a) 
H{r,r') = H(r',r) (196) 

H(r,r') « H(oo,r') + ( 19c ) 

as |r| — > oo, where H(oo,r') is a finite function of r'. 

By integrating eq.(19a) over r with measure J d 2 r and using (19c), we find a 
fourth property: 

J H(\,\')p(\)d\ = -^ (19d) 

For actual calculations, only the one-dimensional restriction H(X, X') is required. 
For A' ~ A, H is dominated by the logarithmic singularity of the free 2d Green's 
function: 

H(X, X') w — In |A - A' | + regular (A ; w A) (19e) 

27T 

Equation (17) is our central tool for evaluating physical quantities using the 
conjugate- field formalism. We shall refer to the last, Wick-expanded factor as the 
'Feynman diagrams', since they go beyond the semiclassical determinant factor* . 

* In the sense that they involve the interaction piece of the quantum action. 
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The only remaining divergences in eq.(17) are UV ones, and they appear in 
two places: the e~ N factor in Z c i ass i C al ( see eq.(5a)), and the determinant factor^ . 
But these two divergences cancel, as we now show. Using eq.(13) and formally 
expanding in powers of the free massless propagator, we find: 



-2-kN J d\p(\)(-^) xx + UV finite 



This is ill defined. But if we use the UV regularization eq.(2a), the S(X — X') term 
in K (r, r') is smeared, and (20) becomes 

det K —I 

ln L . o9 J 1 = Nine + UV finite (20a) 
L det(— a 2 )- 1 

and hence the singular e-dependence drops out of eq.(17), as claimed. 

The remainder of this paper deals mostly with properties of the functions p(X) 
and H(X,X'), especially in the double scaling limit, and with the summation of 
the Feynman graphs in eq.(17). The evaluation of det if, which is done using 
heat-kernel methods, will be reported on elsewhere. 

To get an idea what the expressions for Feynman graphs look like, we record 
the lowest-order terms in the Wick expansion: 

ln( e - Si W) = ln([ J p(X)dXe iq ^] N exp(y J p(X)dXq(X) 2 )} 

= ^{f dr[(q(r f) - [ dr> <g(r') 2 )] 2 (21) 



-2 || drdr'( q (T) q (r')} 2 + ...} 
In eq.(21), the leftmost equality is exact, and we have introduced a useful new 



f Superficially, it appears that the normal-ordering contributions to the Feynman-graph factor 
are also divergent, due to (19e); but these normal-ordering divergences manifestly cancel to 
all orders, as we shall show below. 
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variable r(A): 

A 

r(A) = / p(n)dpL (21a) 



-00 



which ranges from to 1. We will sometimes denote the arguments of a function 
of A, by r instead* . 

The terms displayed in the expansion on RHS of (21), are those involving only 
two propagators. Using eq.(18), we rewrite these Feynman-graph contributions to 
the free energy thus: 

ln(e-^}> =2vr 2 iV{ J dr[H(r,r) - J dr' H(r' ,r')] 2 + 

r f ' ' ' ( 21& ) 

-2 / / dTdT'H(T,r') 2 + 0(H 3 )} 



The first term in the curly brackets is the normal-ordering contribution to this 
order. H(r, r) is logarithmically divergent, but this divergence is r-independent, 
so H{t, t)—J dr' ' H(t' ', t') is UV finite^ . The other two-propagator graph in (21b) 
is manifestly finite. Note that the Feynman-graph contributions to the free energy 
are not local. Indeed, we started from a nonlocal action for the A field, so this is 
to be expected - the cluster expansion does not hold. 

In order to evaluate Feynman graphs, we must solve for the propagator H . 
Before doing so, however, let us show how to sum the normal-ordering contributions 
to all orders. This is easy to do: from the leftmost equality in (21), we obtain 

( e -Sito) = ^J rfre xp(- W (r)) : e iq ^ :} N exp[^ J dr : g(r) 2 :]) (22) 

where w is a finite function, defined as follows: 

w(t) = wi(t) - J dr'wi{T') (22a) 



* In deriving the expansion (21), we made use of the fact that / dTq(r) = 0, by eqs.(14),(21a). 
f This can be shown rigorously by making use again of the UV regularization procedure, eq. 
(2a). 
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wi(t(\)) = lim (2nH(X, X') - In |A - A'|) (226) 

A' — > X 

wi(t) is finite by virtue of eq.(19e). 

We now turn to the task of computing if (A, A'). As was the case with the 
2d equation of motion (7), the two-dimensional differential equation (19a) can be 
converted into a one-dimensiona/integro-differential equation, by treating the RHS 
as a source term. This new equation is 

^H(\ A') = ^T^7 + 2 ^A (p(A)tf (A, A')) (23) 

where the principal value is understood in the first term, and the subscript to the 
Hilbert transform indicates that the transform acts on A, with A' held fixed. When 
the conditions (19b-c) are imposed, eq.(23) has a unique solution. 

It is possible to derive from (23) a Id differential equation for H, similarly to 
the procedure that lead from eqs.(9) to eqs.(10). Namely, eq.(23) is differentiated 
w.r.t. A, then used again, and the properties of the Hilbert transform (listed in 
the appendix) are used. In addition, eq.(9b) for p is used. The result of these 
manipulations is"'": 



{ + 4ir 2 N 2 }H(\, V) = - nN5(r - r') + ^jQ(A|A') 

1 d , p(X') , 



(24) 



27rp 2 (A) <9r /V A- A' 

where principle value is again understood in the last term. Here Q(\\\') is a 
polynomial in A, with coefficients that are one-sided moments of H with measure 
jdr: 

, i i, N V'(X) - V'(X') o f . , , V'(X) - V'U) TT/ xA , v 
Q(A|V) = — ^ _ A/ + 2N J PWn K) x _^ W H(p, X') (24a) 

and these moments are to be determined self-consistently, just as for the moments 
m n of p, which entered the Schrodinger eq.(lOa) through the polynomial -P(A). 



| We use the r variable, defined in (21a), with r' = t(A'), t = r(A). 
§ That is, moments w.r.t. one of the two variables. 
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Denote these one-sided moments as follows: 

h n (fi) = J X n p(X)dXH(X,fi) (246) 

By (19d) we have, 

W = (24c) 
Restricting again to the quartic potential, we obtain: 

Q(X\X') = ^-AgN{[X' 2 + AirNh 2 (X')} + X[X' + 47riWii(A')]} (25) 

The LHS and first term on RHS of (24) constitute the Green's equation for the 
correlator of a quantum-mechanical harmonic oscillator, but the other, nonlocal 
source terms on the RHS spoil this simple picture. As with the eigenvalue density, 
the Id differential equation is somewhat weaker than the Id integro-differential 
equation: the latter, however, is equivalent to the 2d Green's equation (19a), once 
the boundary condition (19c) is imposed. 

The equation (24) is linear, and so can be readily solved in terms of Q (which 
however is itself unknown). The general solution is: 



if (A, A') =A(X') cos 2ixNt + B(X') sin 2nNr - -6{t - r) sin 2nN(r - r') 

2 



(26) 



+ Air 2 N 



o 

where A,B are free functions and 9 is the step function. The symmetry of if (A, A') 
determines A and B up to three real constants; these constants, as well as Q(A|A'), 
can be determined by making use of eq.(23). 

The formalism for evaluating the classical solution(s), propagator and quan- 
tum corrections is rather complicated for finite N . Fortunately, however, massive 
simplifications occur in the expansion, and the formalism simplifies even fur- 
ther when the couplings approach criticality and the d.s.l. (double scaling limit) 
is taken. We next turn to a discussion of this limit. 
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4. The Double Scaling Limit. 

In this section, we will describe the procedure for performing the double scaling 
limit (d.s.l.) in our conjugate- field formalism. We leave out many details, to be 
included in a forthcoming publication. 

Let us specialize to the quartic potential, eq.(lOd), and therefore to the k = 2 
critical model (pure gravity). The plan of the section is as follows. In part 4. a, 
the nature of the planar limit and the d.s.l. for the model is reviewed. Then the 
integro-differential and nonlinear Schrodinger equations for a classical solution, 
p, are used to find the string-perturbative expansions for the moments m n of p 
(defined in (10f)). These expansions are unique, and it is seen they are unique for 
any critical potential, provided attention is restricted to classical solutions that are 
single-band in the planar limit. In addition, the WKB approximation for p in the 
region exterior to the Dyson sea, is found. 

In part 4.b we discuss the perturbative corrections to p inside the sea, and 
study the details of p(A) in the transition regions at the edges of the sea. In 
part 4.c we describe how the same techniques, when applied to the propagator H, 
yield the perturbative expansion for H(\, A') in various regions of the (A, A') plane. 
The formulae of section 3 for quantum corrections are seen to greatly simplify in 
the d.s.l, allowing the perturbative series for physical quantities (specific heat, 
etc.) to be found. In particular, the normal-ordered Feynman diagram expansion 
terminates after a small number of terms. 

Since the critical quartic potential is unbounded from below, all our results up 
to this point were obtained by continuing from the well defined g > regime, to 
g Ri g c < 0. This procedure is satisfactory only for string perturbation theory. In 
part 4.d, we look at the Schrodinger potential V\ directly for g g c ; V\ is seen 
to acquire a small 'second sea' in the transition region. This sea is finite in shape 
when V\ and A are appropriately double-scaled. 

Our formalism is ill defined in this critical coupling regime, like pure gravity 
itself, since there is no normalizable p(A). Nevertheless, we investigate the behavior 
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of p in the new transition region, and find that the second sea has a population 
suppressed by the expected exp(— const ^), where k is the string coupling. The 
constant in the exponential agrees with that obtained in other approaches. Thus, p 
can be thought of as an instanton. As such, it has two unusual aspects. Firstly, the 
tunneling factor occurs in the field configuration as well as in the classical action. 
Secondly, a single conjugate-field configuration seems to describe both the non- 
tunneling and tunneling eigenvalue configurations* . This configuration, however, 
should be viewed only as a warm-up exercise for non-perturbative calculations 
in our formalism. The meaningful calculations are to be done for a potential V 
bounded from below, and would thus most likely not apply to any model with 
k = 2 behavior [29]. 

4. a. WKB and Perturbation Theory for p. 

From eqs.(10), we find the Schrodinger potential for the case of quartic V^: 

Vi(A) = Ag 2 (\ 2 - a 2 )(A 2 - b 2 ) 2 - AgSm 2 + ^(1 + 12,A 2 ) (27a) 

where: 

a 2 = i(vOT%-l) (276) 

b 2 = -— - -a 2 (27c) 
Ag 2 K J 

and 5rri2 is the deviation of the second moment mi from its planar-limit value: 

^^2 + ^-^(1 + 48,) (27d) 
It will shortly be seen that 5m,2 is O(j^) for fixed g, in the planar limit. 



* We thank S. Shcnkcr for discussions concerning this latter point. 

f Our notation, in particular our definition of a(g), differs slightly from that of ref. [46] . 



22 



For g > 0, b 2 is negative, so Vi(A) vanishes only at two points, which are 
A = ±a in the planar limit; these points are the edges of the Dyson sea. We define 
b such that Im b > 0. The k — 2 critical point is at 



g = g c = -1/48 (28a) 

and for g g c , o 2 is positive, and in fact b 2 xs a 2 xs 8 . As explained above, we will 
employ throughout most of section 4 (except part 4.d) the well-defined procedure 
of calculating at g > and then continuing to criticality. The double scaling limit 
for this theory, consists in simultaneously letting N — > oo, g — > g c while holding 
the string coupling fixed [1-8] : 

9 string = K = 1(0 - # c ) _5/4 (286) 



Next, consider our nonlinear Schrodinger equation, (10a), in a region of A 
outside the sea; more precisely, |A| > a, and |A ± a\ are held fixed as N increases. 
The 7r 2 p 2 (nonlinear) term is then exponentially suppressed, by virtue of eq.(9e), 
and in consequence p~ 1//2 approximately satisfies a linear Schrodinger equation: 



1 d 2 

[-jpgx 2 + ( p ~ 1/2 ^ ~ ^ outside sea ) ( 29a ) 



This is easily solved via an asymptotic WKB expansion: (V\ > outside sea) 

A 

p(A) =\VVi(Xj\ exp(-2iV J d^^JT)) {const + 0(1) (2%) 
+ O (e- 27V 1 d ^ v 7 ^) ) } ( outside sea ) 

where Vi(A) is given by eq.(27a)"' : . The corrections on the RHS of eq.(29b) are 
of two kinds: the O(^) corrections constitute the usual, linear- WKB asymptotic 



| It is easy to check that (29b) agrees with eq.(9e). 
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expansion, whereas the exponentially-suppressed corrections are due to the nonlin- 
earities of the exact (10a). Since our V(X) is symmetric, so is V\, and the branch 
we choose for y/V\(p) in the exponent, is as follows: 



sgnyVi(A) = sgnA (29c) 

This choice ensures that p(A) is symmetric, and in addition can be continued to 
an analytic function throughout the complex A plane, except for a cut along the 
Dyson sea (—a, a). 

Next, we take |A| >> I. The exponentially suppressed corrections in (29b) can 
be neglected, and we obtain from eqs.(27a),(29b) an asymptotic expansion for the 
logarithmic derivative of p: 



p'{\) 25m 2 1 l + 12c/A 2 



Np(X) (A 2 -6 2 )VA 2 -a 2 %iV(A2-62) v / A 2_ a 2 

1_ A_ 2A J_ 

+ N [ \ 2 -a 2 + \ 2 -b 2)+ { N 2) 



(29d) 



This expression can now be expanded as a Laurent series in j, and compared 
with the corresponding expansion resulting from eqs.(9) to yield the eigenvalue 
moments, m n . The easiest way to compare the two series term by term, is to 



continue both of them to complex A, with \/A 2 — a 2 defined as described after 
(29c). We then multiply both (9b) and (29d) by A n , n any integer, and integrate 
over d\ along a closed contour with large |A|. For negative n, the coefficients agree 
identically. In addition, all odd moments vanish by symmetry^, whilst for even, 
positive n values we find: 



a a 

J® f \ n d\(\ 2 - b 2 )Va 2 - A 2 + — [ \ n d\ 
TT J y ,W 8ngN J 



8gN5m 2 - (1 + 12#A 2 
(A 2 -6 2 )v / ^ 23 A 2 



—a —a 

a n 1 
+ ^ + d b n + O(j^) (neven) 

(30a) 



§ This symmetry can be shown to hold to all orders in 1/N perturbation theory. 
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with 

For n = and 2, this just reproduces eqs.(9d) and (27d) respectively, up to 
0{l/N 2 ) terms, so we gain no new information. For higher n, eqs.(30) give us 
all the moments in terms of 5m,2, to order 1/JV. How is 5m2 to be determined, 
then? it is clear that do must vanish to this order in 1/N, since at N » 1 the 
support of p(X) is the Dyson sea, or very near it, and the b n term on the RHS of 
(30a) cannot occur for such a distribution, thus do = O(^), which gives us the 
requisite missing information . — 



N5m 2 = 1 - bVb 2 - a 2 + O(^) (30c) 

so 5rri2 is indeed 0(1/ N). Substituting this back into (30a-b) yields all nonvanish- 
ing moments m n , to order 1/N. 

It is straightforward to extend this technique to any desired order in 1/N, and 
to take the d.s.l. limit (28b), as well. We thus see that the perturbative genus 
expansion for p(A), or at least for the set of all its moments, is unique and can be 
easily determined, as claimed. Furthermore, the technique extends to any potential 
V, as long as we have a planar limit of p to expand about. When this limit is 
restricted to have a single band, it is unique, and so the perturbative expansion 
about it will also be unique. 



4.b. Sea Interior and Transition Region. Consider the sea interior, 
namely the region |A| < a with a — \X\ fixed (in either the planar- or the double- 
scaling limits). The nonlinear Schrodinger equation (10a) can be solved via the 

The ^ contribution in (30a), comes from the edges of the sea; the remaining two terms 
come from its bulk. 

* If d 7^ 0, H(p(\)) must have poles at A = ±6, to first order in 1/N. This is actually 
impossible for any p(X) for which the Hilbcrt transform is well-defined, since b is imaginary. 
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WKB approximation in this region, as was done in the exterior region. In the 
interior, however, the linearized WKB is of no use. This is because the planar 
limit we wish to expand about is given by eq.(lOg), and thus the nonlinearity is 
crucial here. 

The correct procedure is as follows. Defining new variables , 



A 

t = N J y/-V!(ji)dfi (31a) 

—a 

p(\) = ^y^V[(Xjf(t)- 2 (316) 
we find the differential equation, 



7T 



f" + f-f = 0(^) (31c) 



In the planar limit, the solutions of (31c) are elliptic functions. 

These planar solutions have two free real parameters; this is just the ambiguity 
discussed at the end of section 2, and is resolved by the integro-differential equation 
and the consistency conditions. Let us see how this works. Since we expect f(t) — > 
1 in the planar limit (by (lOg) and (31b)), we need only consider / ~ 1; then (31c) 
informs us that 



f(t) = l + ecos2t+(^- i cos 4t)e 2 + 0(e 3 )+ O(^) (32a) 



where e is a small unknown oscillation amplitude, and t = t+ip, with ip the constant 
phase of the oscillation. By employing the consistency conditions for tuq = 1 and 



** Recall that V\ < in the sea interior. 
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777.2, it can be shown that 



* = O(^) (326) 



and this turns out to mean that the oscillations can be ignore in the d.s.l. ; thus 
we may use 

/ ~ 1 (32c) 

Next, we study the transition regions at the edges of the sea. By symmetry, it 
suffices to investigate the A ~ a transition region. Since Vi(A) has a first-order 
zero at A = a in the planar limit, we may approximate it by a linear function in 
the transition region, since the width of this region vanishes in the planar limit. 
We rescale A and p as follows: 

A - a = N- 2 / 3 [8ag 2 (b 2 - a 2 ) 2 ]^ (33a) 
P = -N-^[8ag 2 (b 2 -a 2 ) 2 ]^ V (y)- 2 (336) 

7T 

The eq.(lOa) then becomes 

d 2 

+ y + r]- 4 )r]^0 (33c) 



dy 2 

The boundary condition for this differential equation is furnished by the approx- 
imate solution inside the sea, eqs.(31b) and (32c), which become in terms of the 
rescaled variables, 

V(V) ~ (-y)~ 1/4 aty<0, \y\ » 1 (33d) 

As N — > oo, the rescaled equation (33c) becomes exact, and (33d) becomes an 
exact boundary condition in the y — * oo limit. On the exterior side of the transition 
region, the asymptotic behavior is 

r](y) w const y -i/4 e (2/% 3/2 at y > Q> | y | >> l ( 33e ) 
In agreement with the exterior WKB solution, (29b). 
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For fixed coupling g, the width of the transition region is 0(iV -2 / 3 ). When g 
is, instead, continued to the critical point g c in accordance with the d.s.L, we find 
from eqs.(27) 

P - a 2 = 0((g - gc) 1 ' 2 ) (34a) 

so by (33a) and (28b), the width of the transition region is of order (g — g^ 1 / 2 /? 2 / 3 . 
Using eq.(33b), we also find that the normalized eigenvalue population in the tran- 
sition region, is of order 

J pdX = 0(1), (346) 

transition 

either for g fixed, or in the d.s.L. This means that of the original N matrix 
eigenvalues, on the order of one eigenvalue are likely to inhabit the transition 
region. 

4.c. Quantum Corrections in d.s.L In part 4. a, we used both the Id 
differential equation and the integro-differential equation for p(A), to find the per- 
turbative expansion for the moments m n of p. A similar procedure can be employed 
for the propagator H(X, A'), by using the Id Green's equation (24), and the corre- 
sponding integro-differential equation (23). In this case, one finds 1/N expansions 
for the one-sided moments h n (\), defined in (24b). The moments h\ and h<i can 
then be substituted in eq.(25). We find, to the leading approximation for large N, 

W -s>^) (35) 

This, in turn, can be used in eq.(26). The unknown functions A,B are determined 
as explained in section 3. 

The results of this analysis, are as follows*. When A and A' are both interior 

* The details will be presented in a separate publication, as will the computation of the 
determinant factor det K. 
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to the sea, and A 7^ A', we have 

H(X, A') = 0(— ) |A| < a, \X'\ < a (36a) 

The function u>i(A) (eqs.(22)), the regular piece of H(X, X') as A' — > A, is 0(1), and 
in the sea interior it is approximated thus: 

wi(X) ~ lnp(A) + const (|A| < a) (366) 

In other regions of A and A', these functions have different behaviors. For instance, 
when A, A' are both in the same transition region, 

H(X, X') — 0(1) (X^a,X'^a) (36c) 

However, the contribution to Feynman diagrams from vertices in a transition re- 
gion, is still suppressed by a power of 1/N for each such vertex, due to eq.(34b). 
Combining the behavior of p and H in the various regions with eq.(22), we find 
that only the first few Feynman diagrams survive in the d.s.l. We are referring to 
diagrams resulting from contractions among normal-ordered vertices; recall that 
an infinite number of normal-ordering contractions have been summed to obtain 
eq.(22). 

4.d. Instanton Configuration in Direct d.s.l. Rather than continuing 
V(X) to criticality after solving for p and H, it is instructive to attempt taking 
the d.s.l. directly This will demonstrate what nonperturbative effects look like in 
the conjugate- field formalism, although a trustworthy nonperturbative calculation 
requires a potential V(X) bounded from below. 

When we use the quartic potential (lOd) with negative g, b 2 is positive (part 
4. a). As g approaches g c from the origin along the real axis (g ~ g c , g > g c = 
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— 1/48), b approaches a thus: 



b>a , b-a = 0{{g-g c ) 1 ' 2 ) (37) 

Thus by (27a), V\ has two small 'seas', concave regions just outside the main Dyson 
sea, where it is negative. Unlike V(A), V\ is bounded from below. Thus we can 
attempt to solve the nonlinear Schrodinger equation (10a) near criticality, ignoring 
the fact that the solution will not solve the integro-differential equation^. We shall 
continue to use eq.(30c) for 5m,2, in eq.(27a). the justification is that, assuming 
oscillations inside the sea are still suppressed (eq.(32c)), substitution of (31b) in 
(lOf) for n = and n = 2 indeed yields (30c), at least to the approximation needed 
to take the d.s.ft . 

We concentrate on the behavior of Vi(A),p(A) in the new transition region; we 
invoke symmetry again and concentrate on the A ~ a region. To that end, we 
magnify the region via a new rescaling, different from (33): 

\-2V2 = AV3(g - g c ) 1/2 x (38a) 
P = -(9 ~ 9cf /4 ax)- 2 (386) 

71 

where 2v^2 is the value of a at g — g c . The Schrodinger equation in this region 
assumes the form, 

r 1 o d 2 ,64y/8. rxl 1 x2 2^ A 4l , . . 

["48 K 5? + ^ X + ^ X - Vf ' - K S^ ] + C 1C " ° <38C) 

This becomes exact in the d.s.l. The old transition region occurs at x + \[2 = 
0(k 2 ^), and is thus part of the new transition region. We are assuming that the 
string coupling k is small, in order to isolate the leading nonperturbative effect. 



f This is similar to the Marinari-Parisi approach to rendering the even-k models well-defined. 

| This argument, as opposed to the one in 4. a, does not use information about large |A|. Such 
information cannot be trusted, as V is unbounded from below there. 

§ The transition region contributes 0{\/N) to both mo and mi, so to eliminate these contri- 
butions we computed the subtracted TO2 — a 2 m . 



30 



To the right of the old transition region, V\ is positive for 

-V2 < x < 

and p is exponentially suppressed. The extra new minimum of Vi is at x = 1/V2, 
and the new sea surrounding it, where V\ < 0, has width Ax = 0(k 1 / 2 ) and depth 

The resulting d.s.l. solution for p in the new transition region, has the following 
properties. Between x ~ —y/2 and x ~ V\ > and p tunnels in accordance 
with (29b). p is thus suppressed in the small new sea, with the WKB suppression 
factor (up to prefactors which can be found) 

A 2 = exp (_i^ (48) 5/4) (39a) 

which is precisely the tunneling factor appearing in the nonperturbative ambiguity 
for pure gravity, using the various previous approaches [19]* . 

In the new sea, we find that ((x) has the following form: 

z 

C(x ) ~ (1)1/4^-2^ (^-1 _ rf2K l/2 A J e ^ di) (396) 



where z is yet a third rescaled eigenvalue variable, appropriate for the extra sea: 

x- -j= = /t 1 / 2 (48)- 3 / 8 ^ (39c) 

When the d.s.l. limit is taken, k small and z held fixed, the corrections to (39b) are 
higher powers of k. In eq.(39b), d\,d2 are two pure numbers, obtained by matching 
(39b) at large and negative z with the WKB approximation, (29b), which holds 
between the two seas. 



★ In this connection, note that our normalization for k differs from the one usually employed 
in the literature. 
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The configuration given by eqs. (38-39) can be interpretted as an instanton. For 
well-defined models, there exist similar instantons which are true classical solutions 
of the conjugate-field theory. 

In the configuration discussed above, the center z = of the new sea is a local 
minimum of p. p then increases with z for z > 0. When z is sufficiently large, 
one exits the new sea and enters another V\ > region, where the exponentially 
suppressed terms in eq.(29b) dominate for a while; this allows p to continue to 
increase. Eventually the dominant term will again dominate, but if p 2 has by then 
increased to become comparable in magnitude to V\, a nonperturbative solution of 
(10a) is needed. We need not worry about this exterior region, however, since in 
this model p is not a trustworthy configuration there. 



5. Conclusions. 

We have presented a new field-theory formulation of D = matrix models. 
The field is conjugate to the Jevicki-Sakita collective field, i.e. conjugate to the 
density of matrix eigenvalues. The theory is two dimensional, with an eigenvalue 
coordinate and an auxiliary coordinate that can be eliminated from the formal- 
ism. The action is nonlocal, but the equation of motion is local, except for a 
self-consistency condition. There is a unique or discretely labeled classical solution 
for any well-defined potential. The equation for the classical eigenvalue distribu- 
tion is a modified version of the planar integral equation of Bessis, Itzykson and 
Zuber, with an entropy term that smoothes the edges of the Dyson sea and in- 
troduces higher-genus effects already at the classical level. Single-band classical 
solutions are perturbatively unique. The classical distribution also satisfies a non- 
linear Schrodinger equation, with a potential similar to, but different from the one 
appearing in the Marinari-Parisi D — 1 reformulation. 

The classical solutions, and the quantum corrections about them, are systemat- 
ically calculable, and all divergences (UV- and IR-) cancel manifestly to all orders. 
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The normal- ordering graphs can be summed exactly to all orders. In the double 
scaling limit, the formalism simplified drastically. 

The conjugate- field formalism can be used to systematically compute string- 
nonperturbative effects. We demonstrate this for the ill-defined, but simple, case 
of k — 2 realized with a quartic potential. In this case, the classical distribution 
contains two small seas on either side of the main Dyson sea. The population of 
the new seas is exponentially suppressed by the same tunneling factor as appears in 
other approaches. A more complete presentation of the conjugate-field formalism 
will be presented in a forthcoming publication. 



APPENDIX 

In this appendix, we derive the Schrodinger equation (10) from the integro- 
differential equation, eq.(9b), and prove that a classical solution does not have any 
normalization preserving zero modes. 

First, we list a few useful properties of the Hilbert transform (which in our 
normalization is given by eq.(9c)). For any function /(A), 

W(A)) = ^«(/(A)) (Al) 

Also, for any two functions u and v, one easily finds*: 

H{uH{v) + vH{u)) = H{u)H{v) - n 2 uv (A.2) 

Actually, for the present derivation, we only need the degenerate form of this 
identity when v = u: 

2H{uH{u)) = H(u) 2 - vrV (A3) 

But the more general (A.2) is needed to derive the Id differential equation for the 
propagator H(\, A') , eq.(24), from the corresponding integro-differential equation 

* By using the fact that the real and imaginary parts of an analytic function, having all its 
singularities in the lower half of the complex plane, are related by the Hilbert transform. 



33 



(23). We will not go through this latter derivation, but the procedure parallels 
that for p. 

Differentiating eq.(9b) w.r.t. A and using (A.l), we obtain 



Now we use (9b) again to rewrite p' on the RHS of (A. 4), and also use eq.(A.3). 
Upon rearranging terms and using the definition (9c), the nonlinear Schrodinger 
equation, given by eqs.(10), is obtained. 

Next, let Sp be an infinitesimal zero-mode of p which preserves the normaliza- 
tion condition (9d). Thus 

J d\5p(X) = (A5) 
and by varying (9b) we find a linear integro-differential equation for p: 

I(^)' = 2H(M (A6) 

Let us define the function n(r) in the two-dimensional space, as follows: on the 
eigenvalue axis, it is defined as 

n{\) ee S -P (A7) 
P 

and off that axis, it is defined by continuation of (A. 6): 

rj(r)= const + 2N J n(p)p(p) In |r - r'\dp (A8) 
where r' = (p, 0) runs over the eigenvalue axis in the integral. The additive 
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constant in eq.(A.8) is determined by the condition (A. 5), which we rewrite as 

J r](X)p(X)dX = (A9) 

From (A.8) and (A. 9) we find that 

n(r) — > const + Oir—r) as \r\ — > oo (A 10) 

|r| 

But by eq.(A.8), n satisfies the 2c? linear differential equation 

d 2 n{r) = 4irN5(x)p(\)r](r) (All) 

where r = (\,x). This differential equation has a unique solution subject to the 
boundary condition (A. 10). We have thus proven that a classical solution has no 
normalization-preserving zero modes. 
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